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Low-lying optical phonon modes in the filled skutterudite CeRu4Sbi 2 
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The phonon dynamics of filled skutterudite CeRu4Sbi2 have been studied at room temper- 
ature by inelastic neutron scattering. Optical phonons associated with a large vibration of Ce 
atoms are observed at a relatively low energy of E ~ 6 meV, and show anticrossing behavior 
with acoustic phonons. We propose that the origin of the low lattice thermal conductivity in 
filled skutterudites can be attributed to intensive Umklapp scattering originating from low-lying 
optical phonons. By an analysis based on a Born-von Karman force model, the longitudinal 
force constants of the nearest Ce-Sb and Ce-Ru pairs are estimated to be 0.025 mdyn/A, while 
that of the nearest Ru-Sb pair is estimated to be 1.4 mdyn/A, indicating that the Ce atoms 
are bound very weakly to the surrounding rigid RuSb6-octahedron cages. 

KEYWORDS: filled skutterudite, rattling motion, phonon, inelastic neutron scattering 



The rattling phenomenon, that is, a large vibration of 
an atom in an oversized atomic cage, is a key issue in un- 
derstanding exotic physical properties in a family of bulk 
materials that possess cage-like units in their crystalline 
structure. Filled skutterudites, RM4X12 (R = rare-earth; 
M = Fe, Ru or Os; X = P, As or Sb), are one of those 
compounds that have large Xi2-icosahedron atomic cages 
filled with rare-earth atoms. 1 4 Despite their well-defined 
crystalline structures, their lattice thermal conductivity 
is remarkably small, being comparable to that of vitreous 
silica. 5 It has been pointed out by structural refinement 
studies that the thermal parameters of the filling rare- 
earth atoms are unusually large, indicating a large vi- 
brational amplitude of the rare-earth atoms. 6 This result 
has led to the speculation that the rattling motion of the 
rare-earth atoms may strongly scatter acoustic phonons, 
which carry most of the heat flow in a crystal, result- 
ing in an anomalous suppression of the lattice thermal 
conductivity. 

The rattling motion can also influence the electronic 
properties via electron-phonon coupling. The filled skut- 
terudites show various interesting features such as heavy 
fermion superconductivity 7 and a metal-insulator tran- 
sition, 8 which could be affected by the rattling. There 
are also theoretical works that indicate a correlation be- 
tween electronic properties and the motion of filled rare- 
earth atoms. 9 Recently, it has been proposed that rat- 
tling may assist the appearance of superconductivity in 
the /3-pyrochlore compound KOS2O6. 10 Charge fluctua- 
tion enhanced by rattling is considered to be a possible 
exotic pairing mechanism responsible for the supercon- 
ductivity. 

Although the importance of rattling phenomena in 
skutterdites is recognized, there have been only a few 
studies of rattling motion. 11 16 Powder neutron scatter- 
ing and heat capacity measurements suggest that the 
rattling motion is an incoherent localized mode, which 
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can be well described by a localized Einstein mode. 11 ' 12 
The energy of the vibrational rattling motion has been 
estimated from the peak in the phonon density of states 
measured by powder neutron scattering to be E = 5 ~ 7 
meV. Ultrasonic measurements of PrOs4Sbi2 suggest 
that a Pr site splits into four off-centered positions, and 
the Pr atoms vibrate among those positions. 13 

To clarify the influence of rattling motion on elec- 
tronic and thermal transport properties, further studies 
are quite important. Inelastic neutron scattering using 
a single crystal is a powerful method of studying rat- 
tling since it can determine both the energy and mo- 
mentum dependences of phonon spectra. In this work, 
we report an inelastic neutron scattering study that uses 
single crystals of CeRu4Sbi2 and demonstrate the in- 
terplay between acoustic and low-lying optical phonon 
modes characterized by the large vibration of rare-earth 
atoms. We discuss the environment of the filled rare- 
earth atoms based on the analysis using the Born-von 
Karman method and a possible mechanism of extremely 
low lattice thermal conductivity in filled skutterudites. 

Single crystals of CeRu4Sbi2 were grown by an Sb flux 
method, as described elsewhere. 17 The size of a single 
crystal is about 2 ~ 3 mm in length with an almost cubi- 
cal shape. To increase the signal intensity in the neutron 
scattering measurements, six single crystals with a total 
volume of ~ 0.2 cc were assembled. 

Neutron scattering measurements were carried out us- 
ing the triple-axis spectrometer, TOPAN, at the JRR- 
3M reactor of JAEA at Tokai. The final neutron energy 
was fixed at Ef = 30.5 meV using a pyrolytic graphite 
monochromator and an analyzer. A pyrolytic graphite fil- 
ter was used to reduce neutrons from higher-order reflec- 
tions. The sequences of the horizontal collimators were 
15'-15'-S-15'-30' or 40'-60'-S-60'-80', where S denotes the 
sample position. The measurements were conducted at 
room temperature. 

For low-energy phonons at below E = 10 meV, the 
observed phonon peaks were fitted using the following 
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Table I. Interatomic force constants obtained from analysis based 
on Born-von Karman model. 
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scattering function convoluted with the resolution func- 
tion, 



S(q,E) 
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where E s , T, fc#, A and q denote the phonon energy, the 
linewidth, the Boltzmann constant, a scaling factor, and 
the wave vector, respectively. In the fitting, we assume 
that E s has a linear dependence on |q| over the range 
of the instrumental resolution. The dynamical structure 
factor F ine i is estimated using the obtained scattering 
function. F ine i is described as 

F inel = Y -|L eM-W d + iG • r d )(Q • e d ), (2) 



where bd is the coherent scattering length of the dth atom 
at rd, rrid is the mass, e~ Wd is the Debye- Waller factor, Q 
is the scattering vector, G is the reciprocal lattice vector, 
and ed is the polarization vector. It is known that Fi ne i is 
related to the energy-integrated scattering function for a 
one-phonon process in the neutron energy-loss mode as 

S(Q,E)dE ex i dF in ei| 2 . (3) 
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To estimate the interatomic force constants, we have 
performed calculations based on the Born-von Karman 
atomic force model. The longitudinal force constants of 
the seventeen closest atomic pairs were chosen as fitting 
parameters (see Table I), and the calculated intensities 
of the phonon spectra as well as the energies were fitted 
to the measured data. The atomic coordinates of the Sb 
atoms at 24g sites in the space group Im3 were assumed 
to be (0,0.34105,0. 15744) 18 for the analysis. 

Figure 1 shows energy spectra observed along q = 
(C? — C 0) below E = 10 meV, which give the transverse 
phonon modes with the propagation vector [110]. The 
solid lines depict the calculated profiles from the fits con- 
voluted with the instrumental resolution. We find that 
the linewidths of these spectra are comparable to the 



instrumental resolution. At £ = 0.15, a well-defined sin- 
gle peak is observed, which can be assigned to a trans- 
verse acoustic (TA) phonon. Around ( = 0.3, on the 
other hand, two peaks are observed with a separation 
of ~ 1.5 meV. The intensity of the lower acoustic mode 
is strong for small q. As ( increases, however, it gradu- 
ally decreases and vanishes near the zone boundary by 
transferring the spectral weight to the higher branch. 
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Fig. 1. Energy spectra of transverse acoustic and optical phonon 
peaks with propagation vector [110]. The solid lines are the re- 
sults of fits convoluted with the instrumental resolution function. 



The dispersion relations of the peaks observed in Fig. 
1 are summarized in Fig. 2(a). As indicated, the TA 
and optical modes show typical anticrossing behavior at 
£ ~ 0.25. The energy of the lower phonon mode increases 
linearly with increasing £ from the T-point like a typical 
TA mode. When the upper phonon mode appears, how- 
ever, the linear relationship breaks down. In contrast, the 
energy of the higher mode increases above £ = 0.2 and 
begins to saturate near the zone boundary. The behav- 
ior of the dynamical structure factor shown in Fig. 2(b) 
is also consistent with the mixing of two modes. In the 
region above £ = 0.2, the dynamical structure factor of 
the lower energy mode is strongly suppressed, while that 
of the higher energy mode is enhanced, indicating that 
they satisfy a sum rule. 

Figure 3 shows typical profiles of the transverse optical 
(TO) modes for energies greater than E = 10 meV ob- 
served at three q-positions with the propagation vector 
along the [110] direction. Two peaks can be recognized 
at E ~ 18 meV and 35 meV, indicating that the optical 
phonon modes over E = 10 meV are classified into two 
bands. The solid lines depict the results of Gaussian fits. 
Despite the instrumental energy resolution of E Tes ~ 6 
meV, the full width at half maximum (FWHM) of the 
phonon profiles at E ~ 18 meV is about AE = 15 meV, 
clearly larger than E res , suggesting strongly that the pro- 
files actually consist of a number of closely lying optical 
branches. 
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Fig. 2. (a) Phonon dispersion curves of transverse acoustic 
and optical phonon modes with propagation vector [110] in 
CeRu4Sbi2- The motion of rare-earth ions in the optical mode 
(guest mode) is indicated by a red arrow in the inset, (b) Dynam- 
ical structure factor as function of (. The solid and dashed lines 
in (a) and (b) depict the results of a fit based on the Born-von 
Karman model. 
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Fig. 3. Energy spectra of transverse optical phonons with prop- 
agation vector along the [110] direction. The solid lines are the 
results of Gaussian fits. 



To interpret the phonon profiles in Figs. 1 and 3, 
we performed analysis based on the Born-von Karman 
model. From fits to the data, we obtained the calculated 
phonon dispersion relations. Those with the propagation 
vector [110] are depicted in Figs. 2(a) and 4(a). In Figs. 
4(b) and 4(c), we separately show the dispersion of the 
longitudinal and transverse modes. The shade scale of 
the colored curves corresponds to the spectral weight of 
the scattering function. We also depict the peak positions 
of the observed phonon profiles using green circles. The 
vertical lines drawn at each data point above E = 10 
meV correspond to the FWHM of each phonon profile. 
Clearly, the calculated lines follow the observed phonon 
peaks quite well. In particular, the anticrossing behav- 
ior of the observed phonon dispersions and dynamical 
structure factor shown in Fig. 2 are successfully repro- 
duced by our Born-von Karman force model analysis. 



We have confirmed that these anticrossing behaviors can 
be reproduced only with small Ce-Sb and Ce-Ru force 
constants. The alternative lowering of any other force 
constant leads to large separation between calculation 
and observation. The results indicate that the optical 
mode in the low-energy region at E ~ 6 meV is a guest 
mode with T u symmetry. As illustrated in the inset of 
Fig. 2(a), the guest mode allows large vibration of the 
rare-earth atoms. Since the flat optical band produces a 
large phonon density of states in the energy spectrum, 
the guest mode can be the vibrational mode observed 
by inelastic powder neutron scattering measurements at 
E = 5 ~ 7 meV, 11 ' 12 which was identified as being 
the rattling motion. Table I summarizes the obtained 
force constants. The evaluated force constants for Ce- 
Sb and Ce-Ru pairs are very small, 0.025 mdyn/A, re- 
flecting the low phonon energy, and supporting the pic- 
ture that loosely bound Ce ions rattle within the cage. 
However, the anticrossing behavior indicates clearly that 
the guest atoms interact with host lattices, although the 
force constants are very small. This indicates that the 
vibrations of the guest atoms can propagate to the next 
guest atoms through host lattices. The dispersive phonon 
modes indicate that the guest modes at E ~ 6 meV can 
be interpreted as coherent optical phonon modes rather 
than a local incoherent Einstein mode. In the interme- 
diate energy group, the optical phonon energy mainly 
depends on Sb-Sb force constants. These force constants 
have medium values. In the highest energy group, the 
optical phonon energy is determined by the large force 
constants of Ru-Sb pairs, 1.4 mdyn/A, indicating that 
the RuSb6 octahedron is quite rigid. 
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Fig. 4. Phonon dispersion curves along [110] direction in 
CeRu4Sbi2- (a) Results of fit based on the Born-von Karman 
model, (b) Longitudinal and (c) transverse phonon dispersion 
curves. The solid circles depict the results of measurements. The 
vertical lines drawn at each data point above E = 10 meV depict 
the FWHM of each peak. The contour maps show the calculated 
scattering function. 



It has been generally believed that the extremely low 
lattice thermal conductivity in filled skutterudites is re- 
lated to the incoherent rattling motions. Rattled rare- 
earth atoms were considered to act as scattering cen- 
ters for acoustic phonons. This scattering process be- 
comes a dominant mechanism when the rattlers are dis- 
ordered and/or vibrate incoherently. However, the Ce 
atoms in CeRu4Sbi2 are ordered and fully occupy the 



4 J. Phys. Soc. Jpn. 



Letter 



Author Name 



Sbi2-icosahedron atomic cages. Furthermore, our results 
suggest that the guest mode is the coherent optical 
phonon branch, and consequently the Ce atoms cannot 
function as scattering centers. In fact, the lattice thermal 
conductivity of the 100%-filled samples is higher than 
that of partially filled samples where rare-earth atoms 
can scatter acoustic phonons due to disorder. 19 Feldman 
et at. have also pointed out that the scattering mecha- 
nism is not valid for filled skutterudites. 20 Clearly, an- 
other mechanism is needed to explain the low lattice 
thermal conductivity of 100%-filled skutterudites. 

One of the most plausible scattering mechanisms in 
filled skutterudites is their unique Umklapp processes. 
Generally, the phonon modes that contribute to the 
Umklapp processes are restricted to a narrow q-range. 
Figure 5(a) illustrates a conceptual diagram of phonon 
dispersion relations in a conventional system. The gray 
area depicts a trace of the parallel translation of the orig- 
inal acoustic branch by shifting its origin (O') along the 
dispersion relation from the zone center towards the zone 
boundary. The point q' depicts the intersection between 
the optical and shifted acoustic branches, that satisfy the 
conservation law. Clearly, the transition of the "acoustic 
+ acoustic — > optical phonons" through the Umklapp 
process is limited to the short thick line. In contrast, 
for filled skutterudites, the flat optical phonon branch 
lies in a low-energy region below the top of the acoustic 
phonon branch (Fig. 5(b)). Consequently, the process of 
the "acoustic + optical — > optical phonon" transition is 
allowed in a wide area, which can be defined by moving 
the point O' along the low- lying optical phonon branch. 
The optical phonons created through this process are 
distributed along the long q' thick line, suggesting that 
the Umklapp process occurs more frequently than in Fig. 
5(a). It should be noted that this process requires a spe- 
cific condition: namely, the energy of the guest modes 
(^guest) has to be larger than the gap energy between 
the acoustic and upper- lying optical phonon modes (A). 
In fact, A for TA phonons in CeRu4Sbi2 is ~ 4 meV, 
which is smaller than E guest ~ 6 meV. As shown in 
Fig. 4(b), this condition E guest > A is also satisfied 
for the longitudinal acoustic mode. At the same time, 
a small value of £g Ue st is a ls° essential since the number 
of phonons contributing to the process increases with 
decreasing phonon energy. Since these conditions are all 
satisfied in the present system, we propose that Umklapp 
scattering can be one of important processes in the sup- 
pression of thermal conductivity in filled skutterudites. 

In conclusion, we have clarified that Ce atoms are 
weakly bound to surrounding atoms, which ensures their 
large amplitude of vibration. The observed phonons at 
E ~ 6 meV are identified as the optical phonons of a 
guest mode. The results suggest that the remarkably low 
lattice thermal conductivity in filled skutterudites can- 
not be due to an Einstein oscillation. As one possibility, 
we propose that it can be attributed to intensive Umk- 
lapp scattering caused by their unique phonon dispersion 
relations. 
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Fig. 5. Conceptual diagram of phonon dispersion relations, (a) 
Simple system without guest modes, (b) Filled skutterudites with 
guest modes lying within acoustic phonon branches. The gray ar- 
eas indicate the range where a phonon can be created beyond the 
Brillouin zone via (a) acoustic-acoustic and (b) acoustic-optical 
phonon scattering. The thick segments of the flat lines indicate 
the q-range of the resultant phonons. A indicates the gap energy 
between an acoustic phonon and an upper-lying optical mode. 
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